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ABSTRACT 

In the attempt to alleviate the difficulties created by their early formation, we study a 
model in which supermassive black holes (SMBHs) can grow by the combined action 
of gas accretion on heavy seeds and mergers of both heavy (m^ = IO^Mq) and light 
(mf = IO^Mq) seeds. The former result from the direct collapse of gas in T!} ^ 
1.5 X 10"* K, H2-free halos; the latter are the endproduct of a standard H2-based star 
formation process. The H2-free condition is attained by exposing halos to a strong 
( J21 ~ 10'^) Lyman- Werner UV background produced by both accreting BHs and stars, 
thus establishing a self-regulated growth regime. We find that this condition is met 
already at z ~ 18 in the highly biased regions in which quasars are born. The key 
parameter allowing the formation of SMBHs by z = 6 — 7 is the fraction of halos 
that can form heavy seeds: the minimum requirement is that /heavy ^ 0.001; SMBH 
as large as 2 x 10^^ Mq can be obtained when /heavy approaches unity. Independently 
of /heavy, the model produces a high-z stellar bulge-black hole mass relation which 
is steeper than the local one, implying that SMBHs formed before their bulge was in 
place. The formation of heavy seeds, allowed by the Lyman- Werner radiative feedback 
in the quasar- forming environment, is crucial to achieve a fast growth of the SMBH 
by merger events in the early phases of its evolution, i.e. z ^ 7. The UV photon 
production is largely dominated by stars in galaxies, i.e. black hole accretion radiation 
is sub-dominant. Interestingly, we find that the final mass of light BHs and of the 
SMBH in the quasar is roughly equal by z = 6; by the same time only 19% of the 
initial baryon content has been converted into stars. The SMBH growth is dominated 
at all epochs z > 7.2 by mergers (exceeding accretion by a factor 2 — 50); at later 
times accretion becomes by far the most important growth channel. We finally discuss 
possible shortcomings of the model. 
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1 INTRODUCTION 

Ob servations of bri s ht quasars at h ig h redshift 
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and references therein) pose serious challenges to our under- 
standing of their central engine. In fact, their large bolomet- 
ric luminosities, ~ 10*® erg s~^, and hard emission spectra 
are inconsistent with a stellar nature of their energy source. 
Hence, the current paradigm assumes that the central engine 
is powered by baryonic accretion onto a central supermassive 
black hole (SMBH). However, even this proposal is not im- 
mune from additional questions, the most urgent of which 
concerns the growth time of the SMBH. Assuming that a 
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fraction (1 — e) of the matter is accreted at the Eddington 
rate, the growth rate of the SMBH can be written as 



dlnm 
dt 



1 - e 1 

e tE 



(1) 



where Ie = 4:nG^imp/aeC — 0.45 Gyr is the usual Eddington 
time. In order to achieve a SMBH of mass m{t) at a cosmic 
time t{z) corresponding to redshift z it is then necessary to 
start from a BH seed of mass 



mo = m(t) exp 



1 - e t{z) 



tE 



(2) 



From the above expression it is clear that assembling the 
SMBH mass (m = 2 x 10^ Mq) recently deduced for the 
most distant quasar ULAS J1120-I-0641 at z = 7.085 
ijMortlock et al.1 I2OIII I when t{z) = 0.77 Gyr, requires 
lnmo/M0 > 21.4 - 1.71(1 — e)/e. For the usually assumed 
value of e = 0.1, this translates into mo > 4OOM0. Such 
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value is uncomfortably large when compared to the most re- 
cent estimates of the mass of first stars, which tend to con- 
verge towards values well below lOOMm (iGreif et al.|[201ll : 
lOmukai et al.i | 2010,: Bromm fc Yoshidall201ll ). Some authors 
i Willott et all I2OI0I ) investigated the possibility that even 
these small seeds (< lOOM©) can produce large SMBH 
meisses (~ IQ^^Mq), provided they undergo several major 
mergers and spend some time accreting at super Eddington 
rates; current observations support the fact that occasional 
events of super Eddington growth can occur at 2 ~ 6, and 
hence, under those conditions, low mass seeds do not pose 
a serious threat to explain such large SMBH masses. Al- 
ternatively, one has to assume that SMBH growth proceeds 
through the merging of heavier seeds; this possibility is ex- 
plored in the present work. The main questions then shift 
to (i) the formation process of such intermediate mass black 
holes, and (ii) their number density evolution through time. 

Recent theoretical developments in the study of gas col- 
lapse inside dark matter halos propose new mechanisms in 
which heavy black hole seeds of ~ IO^Mq can form rapidly 
(At « 1 Myr) , and hence lend support the alternative sce- 
nario above. In brief, gas inside halos with virial tempera- 
ture Tvir ~ lO'' K cools almost isothermally provided that it 
remai ns mol ecule- free (in practice , H2-f ree). lOh fc HaimanI 
(|2002l ) , and iLodato fc NataraianI (|2006l ) have convincingly 
shown that if H2 formation is inhibited, a primordial gas 
disk is stable to fragmentation and a single massive object 
is formed. Hence, a key point for the mechanism to work 
is that the collapsing halo is exposed to a sufficiently high 
Lyman- Werner (10.2 < hv/eV < 13.6) soft UV intensity to 
photo-dissociate the H2 (or the catalyzer H~) via the two- 
step Solomon process. However, while the UV radiation field 
favors the formation of heavy BH seeds, it inhibits for the 
same physical process the formation of stars in minihalos, 
i.e. halos with viral temperature < 10* K, which rely on the 
presence of molecular hydrogen to enable star formation. 

The UV background is produced by both massive stars 
populating early galaxies and miniquasarfl As the Lyman- 
Werner radiation intensity rises following the formation of 
first stars and black holes, it is expected that star for- 
mation is suppressed in minihalos and continues only in 
T^ir > 10* K halos not fulfilling the requirements for the 
formation of heavy BH seeds. This process comes to an end 
when heavy element enrichment of the halos makes gas frag- 
ment ation unavoidable CS chnci dcr ct al. . 2002 ; Omukai et al.l 
I2OO5I : iDopcke et al.ll201ll ). The interplay between UV back- 
ground effects and the outcome of the collapse of the earli- 
est s tructures has been discussed extensively in the litera- 
ture jRicotti e t al."200 2l:lKuhlen fc Madaull2005l:lAbel etakl 



I2OO7I . l2009inYoshida ct al.ll2007l : Ijeon et al 



2011) and can 



hardly be overlooked. The approach in these papers is com- 
plementary to ours. In fact, these studies, mostly for reasons 
related to the limitations set by the limited dynamical range 
of numerical simulations, have looked in detail at the feed- 
back effects due to the presence of a miniquasar onto the 
star formation in the host (mini-)galaxy. Our approach in- 
stead, being based on a statistical large-scale scheme (the 
merger tree), allow us to study the growth of the SMBH in 



^ We gcncrically refer with this term to lower mass, higher red- 
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the presence of the collective UV radiation field produced by 
the galaxies in the QSO environment. Its novelty consists in 
the fact that a fully self-consistent approach to this prob- 
lem implementing the key processes discussed above has not 
been attempted so far. 

In this paper we intend to specifically address the ques- 
tion of the relative abundances of stellar, light black hole 
seeds and the heavier ones originating from head-start di- 
rect gas collapse including the radiative feedback effects pro- 
duced by the UV background. Followin g t he canonical proce- 
dure in the field JVolonteri et al ]|2003al lbl: [Tanaka fc HaimanI 
I2OO9I : iNataraian fc VolonterilbOlil v" we will attack the prob- 
lem through a Monte Ca rlo technique based on the extended 
Press- Schechter Ansatz (|Parkinson et al.ll2008l ) to simulate 
the mass distribution and merger history of dark matter ha- 
los at various redshifts in a statistical way. We set the initial 
conditions for baryons at redshift z — 20 and track their 
evolution to z = 6 coupling analytical prescriptions to the 
dark matter merger tree. 

We calibrate our models to reproduce the observed 
properties of 2 ~ 6 quasars in terms of their final 
SMBH mass and we compare them with the empirical lo- 
cal s tellar bulge - black hole mass relation rribuia e/mBH — 
10^ (iMarconi fc HuntI l2003l : iHaring fc Rb3 12004 see also 
IValiante et al.ll201ll for a thorough discussion). We next pre- 
dict the evolution of the two populations of seed black holes. 

Throughout we work in a flat cold dark matter model 
with a cosmological constant (ACDM) cosmology with the 
cosmological parameters given by the current WMAP7 
(|LarsonlT201ll ') best-fit values; fim = f^DM + ^It = 0.27, 



f^A = 1 — fim, ^bh 

dus/dhik = 0. 



0.02249, h = 0.704, as = 0.8 and 



2 MODEL 

In this Section we present the basic features of our model. 
The evolution and dynamics of the dark matter field is 
treated in a statistical way by the merger history of a certain 
dark matter halo of mass M which collapsed at redshift 2. 
Each merger event between DM halos is characterized by the 
dynamical friction timescale, tmerge, that depends, to a first 
approximation, only on the ratio x = M> /M< of the collid- 
ing hal o masses. To co mpute tmcrgc we use Chandrasekhar 
theory (|Mo et al.ll2010l ') which states that 

77a; 



H{z)U 



0.234- 



(3) 



ln(l -ha;2)' 

where 77 is a circularity parameter encoding the eccentricity 
of the orbit decay. We define an event with tmerge < (At) 
{tmerge > (At)), with (At) being the mean time resolution of 
the computation discussed below, as a major (minor) merger 
event. If the gas contained in the dark matter halo potential 
wells can cool efficiently (i.e. the timescale of the relevant 
cooling processes tcooi is smaller than the Hubble time tn), 
stars and black holes of various masses can form; the dark 
matter merger dynamics subsequently drives the assembly 
of progressively heavier objects. 

2.1 Merger tree 

The merger tree is a computational realization of the ex- 
tended Press-Schechter (P-S) Ansatz, which gives an approx- 
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imate expression for the mass distribution f unction of col- 
lapse d dark matter objects in the universe (see lLacev fc CoQ 
1 19931 for a more detailed treatment). Take a halo of mass Afi 
at redshift zi: according to P-S, the mass distribution of its 
progenitors at some higher redshift Z2 > zi is given by the 
following expression 



dN 
dM2 



(^2 



ZI 



^1 f ^ ^ 



d In 1^12 



din Ah 



where 



1^12 = 



5c{z2) — 5c{zi) 



(4) 



(5) 



and S{M) is the variance of the cosmic density field 
smoothed on a mass scale M with a sharp k filter. Assuming 
a spherical collapse model we have 



fps{y) = W-!^exp — - 

V TT \ ^ 



(6) 



In order to obtain the mass distribution dN/dN'h from a 
numerical algorithm, one needs to specify a mass resolution 
Mres; to keep computational times reasonable, all details 
below Mres are ignored. It is useful to define two additional 
quantities, that are the mean number of progenitors P in 
the mass interval Mres < M2 < Mi/2, defined as 



(7) 



and the fraction of mass F of the final object in progenitors 
below the resolution limit, defined as 



F = 



dN M2 
dM2 Ml 



dM2 



(8) 



To split a halo of a certain mass and re dshift in its pro- 
genitors we follow [Parkinson et al.l (|2008l ) . Such procedure 
starts from the parent halo (Mi,Zi); we then pick a red- 
shift step Az\ such that P ^ 1 in order to ensure that the 
halo has at most two progenitors (binary mergers) at red- 
shift 21 + Azi. Given a random number R in the interval 
Q < R < 1, ii R > P the halo is not split and its mass is 
reduced to Mi(l — F) (masses below Mres are considered as 
accreted gas from the surroundings during A^i). If instead 
R ^ P the halo is split in two and a random value M2 in 
the range Mres < M2 < Mi/2 is generated from the distri- 
bution (Q; the two progenitors are assigned masses M2 and 
A/i(l - F) - M2. 

Albeit this procedure seems to work at a first sight, 
there is a shortcoming: by definition, when one draws ran- 
domly a mass M2 < Mi/2, the mass of the other fragment is 
automatically chosen as Mi(l — F) — M2, i.e. the two values 
are equally probable, in contrast with the probability distri- 
bution predicted by eq. (g]). Moreover, such algorithm uses 
equation Q only for M2 < Mi /2 and ignores its predictions 
for M2 ^ Ml/2: this results in a halo bias towards the low 
mass end. The fact that this procedure tends to underes- 
timate the mass of the most massive prog enitors is a well 
known problem; to overcome this problem [Parkinson et al] 
(|2008D suggest the following modification of the distribution 
function: 



dN 
dM2 



dMa Wi ' CTi 



(9) 



The dependence on the first argument allows the distribu- 
tion of fragments to be modified, while the dependence on 
5i/ai allows the splitting rate to change; only modifications 
at first order in InG are retained yielding 



G 



02 &\ 

(Ti ' (Ti 



Go 



CTl 



(10) 



The values Go, 71, 72 were chosen so to match the code out- 
put to that of the Millennium Simulatiorl3, which is one 
among t he largest cosm ological simulations available to date 
(see also ISpringelllioOsI ). According to this analysis we set 
Go = 0.57, 71 = 0.38, 72 = -0.01. 

Beside the above biasing problem, merger trees have an- 
other important limitation: matter elements with m < Mres 
are actually accreted by the merged halos but are ignored 
by the numerical procedure. These unresolved mass elements 
can account for a non-negligible fraction of the total dark 
matter mass (~ 90% at 2 = 20) of the simulation. We 
claim that all this missing mass can be reintegrated in the 
tree without any consequences for structure formation, since 
this mass is split in very small elements of mass smaller 
than Mres- The reason is that star and BH formation can 
occur only in halos above a certain mass threshold (dis- 
cussed in the remainder of this Section) corresponding to 
M > 10'' M0 > Mres at 2 = 20 and increasing at lower 2. 
Hence, unresolved structures only bring in dark matter and 
gas, but no stars or BHs. 

We now turn to the discussion on how we populate in- 
dividual halos of the tree with black holes and stars and 
we follow their joint evolution. Before we do this though 
we briefly divert to discuss some important features of the 
mechanism leading to the formation of the heavy BH seeds. 



2.2 Heavy black hole seeds 

Sever al authors (iBromm fc Loebl 20031: Loeb fc Rasid 



1994 lOh fc Haiman^ '20021: 'Lodato fc Natar aianI I2OO6I : 



Volonteri fc Rccs 200^; Regan fc Hach nclt 20091) have pro- 



posed that heavy (~ 10^ Mq) black hole seeds can be formed 
from direct gas collapse, perhaps passing through a very 

short intermediate stellar-like phase^ 

However, hydro-simulations bv lBromm fc Loebl (|2003h . 
suggest that heavy seeds resulting from the direct collapse 
of gas clouds can be formed only if the UV background in 
the Lyman- Werner bands is high enough to photodissociate 
molecular hydrogen. In fact, if H2 is present, the gas cools 
rapidly, thus strongly decreasing the gas fragmentation mass 
scale. Such simulations showed that a mass of ~ 10® Mq can 
condense within a radius (~ 1 pc) where they had to stop 
due to resolution limits; such radius is unfortunately still 
much larger than th e Schwarzschild radius of the system. 
iLoeb fc Rasiol \l994 ) arrived at similar conclusions estimat- 
ing the minimum mass the central collapsed object should 
have in order to stabilize the surrounding disk against frag- 
mentation. 

More recently , lOh fc HaimanI (|2002l ) and 
iLodato fc NataraianI (|2006l ). suggested that the for- 
mation of a heavy black hole seed proceeds via the 
gravitational collapse of accretion disks in H2-free halos 
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Figure 1. Schematic representation of the conditional physical 
processes considered 



and investigated the conditions under which these disks 
can sink into the center of the potential well via the 
Toomre instability. We stress that the formation of a disk 
is necessary in order for gravitational collapse to occur, 
becaus e the typical adiaba tic index of the gas, calculated 
e.g. bv lSpaans fc Sil3 (|2006l ) is too low for spherical collapse 
to proceed without fragmentation. IVolonteri fc Reeg (|2005l ) 
tried to overcome the numerical limitations preventing 
the study of the collapse up to the Schwarzschild radius 
through semi-analytical methods. They were then able to 
investigate the small scale physics at the end stage of the 
collapse. The most remarkable result they obtain is that in 
these systems, because of the high mass of the collapsed 
object (3> Mq), by the time nuclear reactions are able to 
proceed gravity is too strong in order for the collapse to 
be halted. When the temperature reaches extreme values 
(~ 10^ K) weak processes with the emission of neutrinos 
occur and this causes catastrophic cooling that leads to 
the formation of a black hole. Hence, it appears that - 
broadly speaking - the formation of heavy seeds from gas 
collapse is physically possible and it is regulated by the 
strength of the UV background governing the cooling ability 
of the gas. Although this statement is hardly conclusive 
and additional study is required to put it on a more firm 
ground, contingent evidence encourages us to adopt it in 
the following analysis. We check that the gas mass available 
in the halo exceeds the heavy BH seed mass = 10^ Mq. 



2.3 Seeding the tree 

We assign to each halo of mass M on the highest level of 
the tree (2 — 20) a baryonic mass 



We also assign this baryonic mass to all the halos in the 
lower levels that have no parent, that is to say all those 
halos that are too small to undergo further splitting. After 
setting this, we need to decide in which of these halos the 
gas can cool efficiently and hence stars or black holes can be 
formed. Gas gravitational collapse inside the host halos can 
occur only if the cooling time is shorter than the Hubble time 
(tcooi < tn)- Radiative cooling is provided by two distinct 
channels: (i) atomic H i lines, or (ii) H2 molecule rotational 
bands; whether one or the other cooling channel is enabled, 
depends on intensity of radiative background in the Lyman- 
Werner (LW) band ({hv) ~ 12.42 eV). Using the standard 
normalization = J21 x 10~^^erg s~^cm~^Hz~^ sr~^, one 
can define a critical value, J21, separating the two weak and 
strong field regimes. These are discussed separately in the 
next; the seeding prescriptions are graphically summarized 
in Fig. [1] 



2.3.1 Weak UV field limit 

If J21 ^ J21 cooling can proceed only in halos whose mass 
exceeds a threshold value Merit , given by 



10 



-, 3/2 



M.„.(^) = 1.5xlO«rt/eL^^200K;Vl + ^. 

(12) 

iTanaka fc HaimanI (|2009l ') suggest to use the value — 
1200 K-Physically, the excitation temperature of the rota- 
tional levels of H2, the main cooling agent at zero metal- 
licity, would give > 600 K; ho wever a more detailed cal- 
culation (e.g. see IMo et aLll2010l ') yields > 1200 K. The 
critical mass is essentially set by the ability of minihalos to 
cool via the self-synthesized H2. In this case, star forma- 
tion proceeds with efficiency /*. It is yet unclear if a single 
massive ( IOOM0) star forms or rather the final configura- 
tion is a normal stellar population. Although uncertainties 
remain (see Introduction), a common opinion is emerging, 
supported by several lines of observational evidence, that 
the first stars were not as massive as previously thought. 
We take a conservative approach assuming that a stellar 
cluster of total mass nistar ~ f*mb is formed according to a 
Salpeter Initial Mass Function (IMF) 

(/)(m) oc m~", (13) 

with a = 2.35. Stars leaving behind BH remnants are those 
in the mass range 40Mq < m < I4OM0 and m > 26OM0 
(in the intermediate range the star ends its life as a pair- 
instability supernova). The lower cutoff rricut for the frag- 
ments mass fixes the normalization constant and so we can 
integrate eq. (|13|l to obtain a value for fBH^ ^.s follows 



Jbh (iricut) 



(j>[m)dm 



(14) 



i Lfri 



(11) 



Integration yields the following values: /^^^(2OM0) — 0.35 
and /^^^(IMq) = 0.006; the cut at 2OM0 is suitable for 
collapsing gas inside metal-free primordial halos; mcut = 
IMq applies to halos which have been polluted by metals 
from previous generations and have a smaller fragmentation 
scale. For simplicity, we assume that every halo which has 
hosted stars at the previous (i.e. higher redshift) merger tree 



step is metal polluted and therefore f^u 



0.006 
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2.3.2 Strong UV field limit 

If J21 J21, molecular hydrogen is photodissociated and 
cooling can proceed via atomic H lines (and metal lines, if 
present). If cooling is effective (i.e. a condition requiring a 
halo virial temperature ^ 1.5 x 10'* K) then one must 
evaluate the halo spin parameter, A, to infer the fate of 
the gas. Because the situation is still not well understood 
from the theoretical point of view, we make the following 
simplifying assumption: if A > Ac then the gas settles into 
an accretion disk that fuels the formation of a heavy black 
hole seed of mass ruheavy via Toomre instability, on a short 
timescale, typically 1 Myr. If instead A < Ac then the disk 
does not form, the gas cloud collapses nearly spherically and 
fragments, again forming stars with efficiency /, . We recall 
that this scenario is mainly supported by the evidence that, 
given the typical valu es of the adiabatic index calculated by 
ISpaans fc Silkl (|2006l l , a spherical gas cloud is very likely to 
fragment; thus, the formation of a disk is needed to allow 
the formation of a massive central object. We can translate 
the spin threshold criterion into the fraction of halos that 
can form heavy se eds, f heavy , via- the log- normal probability 
distribution p(A) (|Warren et al.lll992l ): 



10-" 



p{X)d\ 



(15) 



The value of Ac depe nds on small scale details (see 
iRegan fc HaehneltllioogI for a thorough discussion) such as 
the gas temperature, Tgas, the collapsed baryon fraction, fd, 
and metallicity, which are all rather uncertain; we then keep 
Ac (or alternatively fheavy) as an additional free parameter. 
We stress that the heavy seeds are thought to be formed on 
a very short timescale t form. ~ 1 Myr compared to the mean 
time step used {At) ~ t/f/15 with tn ranging from 1.4 Gyr 
at z = 6 to 0.3 Gyr at z — 20. The details of the formation 
hence hopefully should not affect our results. 



2.4 Ultraviolet background 

How strong could be the LW background in the biased re- 
gion in which SMBHs we observe at 2 « 6 are built? A 
minimal and yet robust estimate can be derived as follows. 
There is now strong observational evidence that at z = 6 
cosmic reionization was essentially complete. Thus, one can 
compute directly the mean intensity of the ionizing (cosmo- 
logical) UV background, J^, required to keep the intergalac- 
tic hydrogen ionized by requiring that at least a photon per 
baryon is produced, n.y/ni, — 1. However, this estimate rep- 
resents a lower limit to as it does not account for recom- 
binatio ns occurring in oy erdense regions. Recent numerical 
works (|Mitra et al.ll20l"ll ) suggests to take n^y/n^ = 10. It is 
then easy to show that 

hpc {pH[z)) / 
rup \ m 



A7vJ+(z) = 



(16) 



j+{z) = o.u{nbh''){i + zf 



(17) 



which yields J^(2 = 6) = 0.97{n.y/nb) ~ 10. In the previous 
equations, hp is the Planck constant, c is the speed of light, 
rrip is the proton mass, and {ph{z)) is the mean density of 
hydrogen contributing 76% of the total gas mass. However 



103' 



103' 




100 



1000 

m 



Figure 2. Stellar continuum spectrum for instantaneous star for- 
mation mode, normalized to a star cluster of mass m = IO^Mq 
with metallicity Z = 0.001, and an IMFoc m~^-^^ with rricut = 
IMq. Labels shown the value of the spectral jump, .J~^ / J~ below 
and above 1 Ryd. 



the intensity of LW radiation, which is not absorbed by hy- 
drogen atoms, can be larger by a factor equal to the inverse 
of the escape fraction of ionizing photons f^g], « 10^"'^ mul- 
tiplied by the jump across 1 Ryd in the intrinsic spectrum 
of a typical stellar population (see Fig. [Sjl , which is a func- 
tion of stellar ag^fl for the IMF used here. As we see from 
the Figure, the jump is in the range 2.6 — 60.2 for a young 
stellar population and increases further with age. We can 
safely conclude that close to reionization, we expect average 
LW intensities as > J21 ~ 10^. We however expect the (lo- 
cal) values of J21 be higher in the star forming (and high 
quasar activity) regions we consider being higher than this 
average value, and a more detailed calculation is needed to 
correctly estimate the local LW background. These details 
are explained in the subsequent paragraphs 2.4.1 and 2.4.2. 

In practice, though, one needs to properly follow the 
UV radiation sources during their evolution to obtain a 
more precise estimate. We model the emission of sources 
in the merger tree through a redshift-dependent specific 
physical emissivity e(i',z) = L^(z)(l -I- z)^/Vbox, where 
Vbox ~ Mh/pm is the comoving volume corresponding to 
a linear perturbation of mass Mh that we use to determine 
the simulated volume, and = J]]^ Li,^i is the total lumi- 
nosity produced by the sources in the LW band. For the 
Mh = lO^^M© halo we are simulating the comoving volume 
is Vboo: = 29Mpc^ 

Lyman-Werner photons interact with hydrogen 



^ Results obtained with S TARBURST99 llLeitherer et al.lll999l: 
I Vazquez fc Leithereill2005l) . software available at 
|http : /7mww. stsci . edu/scienc6/starburst99/ | 
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molecules present in the intergalactic medium and are 
absorbed. Assuming a H2 re lic abundance « 10~^ one typi- 
cally gets (jCiardi et al.ll2000l ) an optical deptlfl TH2 ~ 3. We 
assume that inside Vbox the radiation field is homogeneous. 
This is justified by the fact that the mean free path of 
LW photons is much larger than the physical size of the 
simulated volume. The background can be computed by 
properly accounting from redshifted radiation from higher 
redshift sources, via the standard formula: 



47r 



l + z' 
' 1 + z' 



l + z 
l + z' 



dt' 



dz' 



dz 
(18) 

where Zon is the source turn-on redshift. 

Apart from more exotic and rather unconstrained radi- 
ation sources as dark matter particle annihilations/decays 
iMapeUi et al.ll2006l : rValdes et al.ll2010l '). there are two types 
of sources that contribute to the UV radiation background, 
i.e. stars and accreting black holes. 



24.1 Stars 

UV radiation is copiously produced by massive stars, par- 
ticularly if young and metal-poor. For example a metal free 
population of stars with a IMq — 500Mq Salpeter IMF, pro- 
duces lO"*-*^^ H2 photodissociating photons per baryon into 
stars while on the ZAMS (compared to IQ''-^^''"' in H -lonizmg 
photons). Given that the narrow width of the LW band (2.4 
eV), the stellar spectrum is assumed to be flat within it; 
its time evolution (Fig. [2| has been derived through the 
public stellar population synthesis code STARBURST99 for 
an instantaneous burst of star formation with metallicity 
Z = 10~^ ^ 0.05^0 for the Salpeter IMF with the two dif- 
ferent low mass cut-offs rricut = (1, 20)i\/Q introduced above. 



2.4-<^ Accreting black holes 

In comparison to stars, black hole emission is currently less 
understood and we need to make three simplifying assump- 
tions in order to compute it. These are: (i) each heavy black 
hole that is surrounded by a cloud of gas accretes mass at 
the Eddington rate 



rriBH = fdu 



ty 



1 



(19) 



where ts ~ 0.45 Gyr; the radiation conversion efficiency is 
taken to be e = 0.1 and we introduced a duty cycle fac- 
tor fduty ~ 0.5 — 1 to account for t he fact that gas ac- 
cretion may proceed at a lower ra te (|Marconi et al. I I2OO4 
iBrusa et aLll2009l : [Fiore et al.ll2012l '). We note that the pa- 
rameter fduty contains a degeneracy between the quasar 
lifetime and the accretion rate. For example, two accretion 
events with fduty = 0.5, one proceeding at half the Edding- 
ton rate for time interval At, and the other proceeding at the 
Eddington rate for half of the time interval At will result in 
the same net accreted mass. When the black hole has a mass 



^ During the evolution, as J21 rises, most of the H2 will be cleared 
by photodissociation with a consequent decrease of the opacity; 
we neglect this complication. 

^ This implies that stars form on a timescalc t/orm i^t)- 



rriBH we assume that its luminosity coincides with the Ed- 
dington limit LE{mBH) = 1.3 x lO'*" erg s~ ^(msff/ 100 A/q); 
(ii) light black holes do not accrete gas and hence do not 
co ntribute to luminosit y. This is supported by the results 
bv I Alvarez et al.1 (|2009t )') who have showed that this type of 
BHs are characterized by very low 1O~^^M0 yr~^) accre- 
tion rate partly due to radiative feedback and also because 
they spend most of their lifetime in low-density regio n s. Thi s 
setting is different than that of iTanaka fc HaimanI (|2009l ). 
where they include the possibility that light seeds do ac- 
crete, principally due to the fact that these seeds form in 
the center of their host halo. In our setting however, light 
seeds are scattered in the stellar bulge, where the gas low 
density prevents them from accreting signi ficantly. We sh- 
old also note t hat iMilosavlievic et al.l (|2009l ) arrived at sim- 



ilar results as lAlvarez et al. ( 20091 ): both simulations used 



the Bond i-Hoyle-Lyttleton formula to estimate the accre- 
tion rate. IMilosavlievic et al.l (|2009l ) also suggest that even 
feedback-suppressed systems may undergo short episodes of 
potentially super-Eddington accretion, but in this paper we 
neglect this possibility; (iii) All accreting black holes are 
characterized by a non-thermal (x l/u emission spectrum in 
the band 10eV< /ii^ <10keV; this sets the normalization 



L„(mBH) = 



1 LEirriBH) 
3 In 10 V 



(20) 



The above spectra for the two types of sources are then 
used to calculate the value of J21 through eq. (|18|) . Finally 
we need to fix the chosen fiux threshold value, J|i, for 
H2 photodissociation. This value depends on the spectrum 
shape and w e refer to recen t work to obtain a numeri- 
cal estimate. iBromm fc Loebl (|2003l ') find J|i « 10^ for a 
quasar power-law spectrum and J|i < 10^ for a lO"* K ther- 
mal spectrum that approximates well a stellar-dominated 
background. More recently, the quasa r threshold has bee n 
decreased by an order of magnitude (|Omukai et al.|[2008l ): 
such finding has been confirmed an d refined by a more 
carefu l estimate of the LW opacity bv IWolcott-Green et al.l 
(|20lj ). These authors suggest to use the threshold values 
J21 = 10^ for a stellar background and J|i = 4.3 x 10^ 
for a power-law quasar background. In our simulation, the 
actual UV background intensity and shape represents a 
combination of the stellar and accreting black hole con- 
tributions and therefore in such a mixed situation it is 
not obvious which of the two thresholds applies. Lacking 
a more definite insight, we adopted the prescription that 
H2 photodissociation occurs when one of the two thresh- 
olds is reached, independently of spectral shape of the dom- 
inant component. This choice is not crucial in practice, as 
we will see later on, BH radiation is always largely sub- 
dominant with respect to the stellar one. To summarize, the 
numerical parameters we use are (J21, T's , Ts*, n^f , /*) = 
(1-4.3 X 10^ 1200 K, 1.5 x 10* K, IO^Mq, IO^Mq, 0.1). 



3 SIMULATIONS 

In this Section we will brieffy describe the simulation im- 
plementation. The me rger tree is built very similarly to 
[Parkinson et al.l (|2008l ) and it determines the dark mat- 
ter halo mass function evolution with redshift. We follow 
the tree in the redshift interval z G [6, 201 with a mass 
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Figure 3. Merging timescale in units of the Hubble time l(3j as a 
function of the mass ratio x = M> /M< for different values of the 
circularity parameter 77. Also shown is the value H{z)tmerge = 
1/15 corresponding to our mean timestep in a tree with 2 G [6, 20] . 



resolution of Mres ~ W^AIq; the adopted redshift step is 
Az = 0.7, corresponding to a mean timestep {At) = 40 Myr 
(or (HAt) — 1/15). This choice is motivated by two rea- 
sons: the first is that, since the heavy black hole seeds are 
thought to form on a timescale tform = 1 — 10 Myr, we can 
ignore the uncertain formation process details. The second 
is that the timestep should be larger than the timescale for 
a major merger, which cannot be arbitrary small; taking a 
mean circularity parameter rj = 0.2, from Fig. Owe identify 
a; = 4 as the critical mass ratio that discriminates between 
major and minor mergers. M otivated by the recent discov- 
ery of iMortlock et al.l (|201ll ) we included also some runs 
with 20 timesteps in z G [7.08, 20], a mean time resolution 
(At) = 30 Myr and a critical mass ratio a; = 3. With such 
high mass and time resolution we can computationally af- 
ford to compute the merger history of a large final SMBH 
host halo of Ah = IO^^Mq at z = 6. The choice of z = 20 
as the initial redshift is dictated by computational reasons: 
with this choice, we achieve a mass resolution of lO'''Af0, es- 
sentially equal to the critical mass (at z = 20) discriminating 
between star-forming and "dark" halos, in which the gas is 
not able to cool. Pushing the initial redshift to 2; = 30 — 40, 
would require a further increase the resolution to retain the 
ability to discriminate star forming halos. We do not expect 
our assumption of 2; = 20 as the initial condition to affect 
greatly our results because the timesteps become very small 
at high redshift, and hence the BH seeds do not accrete a 
significant amount of mass before that time. 

We track the evolution of baryons in each halo in the 
following different components: (i) a central BH, built from 
heavy seeds merging and gas accretion (see below; we ne- 
glect direct accretion of stars by the BH); (ii) a bulge com- 
ponent, containing stars, light BHs and gas available for 



BH accretion or further star formation, (iii) the satellites, 
i.e. sub-halos resulting from minor mergers, which contain 
heavy and light seeds, stars and gas. Some discussion on the 
merger prescription is in order at this point. During a major 
merger all the contents of the light halo sink into the center 
of the heavy halo. In particular, heavy holes merge, whereas 
light holes, stars and gas are simply added to the bulge; the 
contents of the heavy halo satellites remain the same. For 
minor mergers instead, the heavy halo remains unchanged, 
whereas the contents of the light halo are added to the heavy 
halo as satellites. When two black holes merge, we always 
assume that the merging timescale is controlled by the dark 
matter dynamical friction timescale, i.e. black holes merge 
only during major mergers. Throughout this work we neglect 
the final stage of the merger which is regulated by the grav- 
itational radiation loss timescale. The latter is much smaller 
than all other relevant times scales for heavy seeds mergers, 
but it is unclear whether the same condition holds also for 
light - heavy seeds mergers, due to the final parsec barrier. 
Negle cting this not yet well understood pr ocess (see for ex- 
ample I Nixon et al. l l201ll : |Preto et al.ll2011^ introduces some 
uncertainties in our model, but its influence on our results, 
whose uncertainty is predominantly of statistical nature, is 
likely to be subdominant. In all the cases considered only 
the heavy BHs in the bulges are allowed to accrete at the 
Eddington rate (with duty cycle fduty) until they run out 
of fuel; on the contrary, no accretion takes place in the pe- 
riphery, and hence BHs in satellites do not emit. Multiple 
merger events in the tree have been modelled as a sequence 
of binary mergers occurring in a random order. 

Finally a mechanism that can cause the merging of light 
and heavy holes has been included: the light holes in the 
bulge experience dynamical friction due to stars and can 
sink onto the central BH (if present). The calculations for 
this sinking timescale, along with the definition of major and 
minor mergers, are analogous to those for the dark matter 
merging timescale tmerge so we do not repeat them here. 
Using these prescriptions we set the initial conditions for 
baryons at z = 20 and track their evolution to 2; = 6 (or 
z — 7.08) to verify if a central black hole mass of ~ W^Mq 
can be put together. 



4 RESULTS 

In this Section we show the results obtained from the model 
described above. Predictions concern mostly the case for a 
SMBH hosted by either a Mh = lO^^M© or Mh = IO^Mq 
halo a.t Zf — 6. However, motiv ated by the recent di scov- 
ery of the most distant quasar (|Mortlock et al.ll201ll ). we 
have also performed some runs ending at Zf = 7.08. The 
first point we want to make is that the chosen values of the 
accretion duty cycle fduty and of the critical spin parame- 
ter Ac (or equivalently, in our model, of the heavy seeding 
fraction f heavy) set the order of magnitude of the final black 
hole mass obtained at z/ = 6, as can be seen in Figure U 
If accretion occurs at the Eddington rate {fduty ~ 1), the 
black hole mass is almost independent of f heavy and, beyond 
the sharp cutoff at fheavy ~ 0.01, it sets to a value in the 
range Mbh = 6 — 7 x 10^° Mq. Thus, extremely large black 
holes, exceeding the current estimates, can be formed under 
these conditions. The more realistic case, fduty = 0.5, yields 
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Figure 4. Final black hole mass at 2j = 6 as a function of the 
heavy seeding fraction, f heavy > for two values of the duty cycle, 
fduty Error bars show Icr dispersion around the mean over 100 
merger tree realizations. 
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Mh = 10^'^ Mq, Mres = 


= lO^M© 


zf = 6 


0.5 


8.5 ±0.2 


11.02 ± 0.08 


1.0 


10.6 d= 0.1 


10.7 ±0.2 


zf = 7.08 


0.5 


7.8 ±0.3 


11.00 ±0.08 


1.0 


10.0 ±0.4 


10.8 ± 0.1 


= 


6, Afh = 10" Mq, Mres 


= 2 X IO'^Mq 


0.5 


9.1 ± 0.3 


12.06 ±0.06 


1.0 


11.6±0.1 


11.8 ±0.09 



Table 1. Outline of the results for f heavy = 0.006; the final black 
hole mass Mgu and the mass of the surrounding baryonic bulge 
^bulge s-rc displayed. There is also a run with A/^ = 10^^ Mq. 
Error bars show Icr dispersion around the mean over 100 merger 
tree realizations. 



a final central BH mass at Zf = 6 in tlie 10 M0 range, 
strongly dependent on the value of f heavy Table [T] contains 
a summary of the final black hole and bulge masses obtained 
for various fduty and Zf, for a fixed value of fheavy ~ 0.01. 
Let us now examine in more detail the preferred model in 
which fduty ~ 0.5 and fheavy = 0.01 by considering the pre- 
dicted evolution of the radiation background, the star/BH 
formation rates, and the masses of the various baryonic 
components. Looking at Figure (5] we see that the Lyman- 




10 IT) 20 

z 



Figure 5. Evolution of the UV Lyman- Werner radiation field 
J2i{z). Plotted are the stellar (red lines) and black hole (black 
lines) contributions for Zf = 6 (solid) and Zf = 7.08 (dotted) 
for the case Mh = W^^M^Ja^ty = 0.5, fheavy = 0.01. Error 
bars show Icr dispersion around the mean over 100 merger tree 
realizations. 



Werner radiation field, and hence H2 photodissociation, is 
dominated by the stellar background. The specific intensity, 
J2i(z), rises very rapidly and up-crosses the critical thresh- 
old J21 = 10'' already at z ~ 18, thus triggering the forma- 
tion of heavy black hole seeds. This is particularly impor- 
tant, as it shows that in the highly biased regions in which 
quasars are born, the UV radiation from galaxies is sufficient 
to allow the formation of the heavy seed required to build 
them. In this model the black hole UV background com- 
ponent does never become supercritical. The fact that an 
early epoch of star formation is nec essary for the formatio n 
of heavy seeds was also suggested bv lBromm fc Loebl (|2003l ). 
We note that even if the average cosmological background 
reaches the critical threshhold value ~ 10"^ only at the epoch 
of reionization, the local one in star forming regions is sigif- 
icantly bigger, and reaches the threshold for heavy seed for- 
mation much earlier. 

The early appearence of heavy seeds can be better ap- 
preciated from Fig. [6l where the various formation rates are 
compared. The light BH seeds formation rate has an initial 
peak, it decreases and finally stabilizes to a much lower and 
relatively constant value, 0.2 — O.SM© yr^^, as a result of the 
progressive metal enrichment that forces the light seed for- 
mation to proceed with a lower fg^^ = 0.006, as discussed 
in Sec. (2] As expected the light BH formation rate closely 
tracks the star formation rate; the heavy BH seeds forma- 
tion rate instead shows considerable fluctuations around a 
lower mean value. Once formed the heavy BH seeds grow 
by accretion and merging with other heavy seeds, according 
to the merging prescriptions introduced above. The peaks in 
the stars and light BH's formation rates are a numerical arti- 
fact related to the way in which we modeled metallicity. We 
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Figure 6. Formation rate M of stars (blue lines), light (red), 
and heavy (black) black hole seeds as a function of redshift 
for Zf = 6 (solid) and Zf = 7.08 (dotted) and a model with 
Mh = lO^'^MQjduty = 0.5, fheavy = 0.01. Error bars show la 
dispersion around the mean over 100 merger tree realizations. 



Figure 7. Redshift evolution of the stellar, black holes, and avail- 
able baryons for = 6 (solid) and Zf = 7.08 (dotted), for the 
model with A-4 = W^^MeJduty = 0.5, fheavy = 0.01. Error 
bars show Icr dispersion around the mean over 100 merger tree 
realizations. 



made the following simplifying assumption: if a halo did not 
experience star formation before, we consider it as metal- 
free and we set 20A'Iq as a IMF lower cut. If instead star 
formation already occurred, then we assume that its metal- 
licity is above the critical one and cut the IMF at IMq. 
This is admittedly the simplest possible hypothesis and this 
leads to the peaks in Fig. 6. The correct solution is probably 
smoother but this effect is only a marginal one affecting the 
very phases of the system evolution. 

The final mass of the two compact populations, light 
BHs and the central SMBH, is roughly equal by z = 6, as 
seen from Fig. [7] where we present also the evolution of the 
redshift evolution of the baryonic components. However, a 
nice feature of the results is that both populations are sub- 
dominant with respect to stars by 2 — 3 orders of magni- 
tude. As we will discuss in more detail in a moment, this 
is broadly consistent with observation of the stellar bulge 
- black hole mass relation for high-z quasars. Baryons are 
not completely converted into stars even though star forma- 
tion has proceeded at very sustained rates (> SOMoyr^^) 
for essentially the entire Hubble time; we find that 19% of 
the initial baryon content has been converted into stars. The 
baryonic content is continuously replenished by infalling gas 
(in our scheme, baryons contained in very small, starless 
halos with mass below the resolution mass, Mres). 

The SMBH mass Mbh at 2: = 6 is set by the combina- 
tion and balance between gas accretion (we do not consider 
accretion of stars) and BH mergers along the tree hierarchy. 
A natural question is then which of these two process pre- 
vails along the build-up history. This question is answered 
by Fig.[8]where we plot the fraction of the final SMBH mass, 
/^Macc/ Mbh , gained by either accretion or merging in each 



tree redshift step. These curves has been obtained by track- 
ing the evolution of the central BH mass along the most 
massive progenitor branch of the merger tree and by defin- 
ing the "merged" mass as AMmer = Mbh{z) — AMacc For 
the accreted mass in a timestep At (or in a redshift step 
Az) we used the usual expression 



AMacc (Az) 
Mbh{z) 



(21) 



The SMBH growth is dominated at all epochs z ^ 7.2 by 
mergers, which dominate the fractional mass gain by a factor 
2 — 50 depending on redshift. Between 2 = 7.2 and 2 = 6, 
the relation is reversed and accretion becomes by far the 
most important growth channel, with mergers progressively 
becoming negligible with time. During this late evolutionary 
stage the final black hole mass grows by a factor of nearly 
three. As already shown, these results depend on the un- 
known fraction of halos that can form heavy seeds, fheavy 
If one poses that the local stellar bulge - black hole mass 
relation holds also at the high redshifts of concern here, we 
can take advantage of the empirical observation al constraint 
MBH/Mtuige ~ 10"^ (scc IVafiaute et aLllioTll ) to obtain a 
more realistic value for fheavy ■ Decreasing fheavy values tend 
to improve the agreement between the model results at 2 = 6 
and the local empirical relation (Fig. |9]): the best fit is ob- 
tained for if duty, fheavy) ~ (0.5,0.006). Independently on 
fheavy, though the relation is steeper than linear, implying 
that SMBHs formed before their bulge was in place. This 
trend is in line with the results of recent numerical simu- 
lations investigating the nature of quasar hosts at 2 = 6 
(|Khandai et al.ll201ll ). 
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Figure 8. Fraction of the final SMBH mass, AMacc/ Mbh , 
gained by either accretion or merging in each tree redshift step 
for the model with zj = 6, = f^uty = 0-5, f heavy = 

0.01. 
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Figure 9. Mbh — J^Ibulge relation for a tree realizations with 
Zf = a (thin points) at each tree level for three difi'erent values 
of fhnavyy ^ indicated. For the best fit value fheavy = 6 X 10^'' 
we show also the case with Zf = 7.08 (thick points). The dotted 
line is the local empirical relation Mgu = 0.001A/(,„jgg. 



5 CONCLUSIONS 

In the attempt to alleviate the difhculties created by the 
observed rapid formation of SMBHs, which are already in 
place only 770 Myr after the Big Bang, we have proposed a 
model in which the hole can grow by a combined action of 
gas accretion and mergers of both heavy (m^ = IQPMq) and 
light (mf = 10^ Mq) seeds. We stress that, in our picture, 
the heavy seeds are formed in the center of their host halos 
and hence are able to accrete gas, whereas the light seeds 
are scattered in low density regions and hence do not accrete 
gas efficiently. The former result from the direct collapse of 
gas m Ti* ^ 1.5 X 10" K, H2-free halos; the latter are the 
endproduct of a standard H2-based star formation process 
which also depends on the enrichment conditions of the gas. 
The molecular-free condition is attained by exposing halos 
to a strong ( J21 ^ 10^) Lyman- Werner UV background pro- 
duced by both accreting BHs and, predominantly as we have 
shown, by stars. Quite noticeably, UV radiation in the biased 
region in which quasars live can in principle establish a self- 
regulated joint evolution (schematically represented in Fig. 
1) of star formation and black hole formation/growth: fluxes 
exceeding the critical value J21 favor the formation of heavy 
BH seeds, but they simultaneously depress star formation in 
small halos, and hence their UV photon production. With 
our model, based on a well-tested merger tree scheme and a 
treatment of the Lyman- Werner field intensity produced by 
both stars in galaxies and black holes, we have been able to 
follow all these processes in detail. 

We have shown that the key parameter allowing the 
formation of SMBHs by 2; = 6 — 7 is the fraction of ha- 
los that can form heavy seeds: the minimum requirement is 
that, fheavy ~ 0.01; SMBH as large as 2 x 1O^°M0 can be 
obtained when fheavy approaches unity. The precise value 
of th is parameter is very har d to derive from first princi- 
ples (|Regan fc Haehneltl Boogl ) . as it depends on unknown 
environmental properties (gas temperature and metallicity, 
collapsed baryon fraction). However, the mere existence of 
SMBHs places the above solid lower limit to this quantity. 
Independently on the value of fheavy, though, we find that at 
high- 2 the stellar bulge - black hole mass relation is steeper 
than the local one, implying that SMBHs formed before 
their bulge was in place. The formation of these heavy seeds 
is then crucial to achieve a fast growth of the SMBH by 
frequent merger events in the early phases of their evolu- 
tion, z ^ 7. Their formation is allowed by the fact that the 
Lyman- Werner radiation field intensity rises very rapidly 
and overshoots the critical threshold J21 in the quasar- 
forming environment already at 2; ~ 18; such photon input is 
almost completely dominated by stars in galaxies. This UV 
fiux, on the other hand, considerably quenches star forma- 
tion in mini-halos, but we find that this effect is negligible 
as far as the UV intensity is concerned, i.e. weak feedback 
strength. 

Interestingly, the final mass of light BHs and of the 
SMBH in the quasar is roughly equal by z = 6; by the 
same time only 19% of the initial baryon content has been 
converted into stars. The SMBH growth is dominated at all 
epochs z > 7.2 by mergers (exceeding accretion by a factor 
2 — 50); beyond that epoch accretion becomes by far the 
most important growth channel. 

In spite of the novel results obtained, the present work 
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should be seen only as a useful guideline for future and 
more in-depth work. A number of additional physical effects 
should be included along with a more detailed description 
of many aspects of the problem. 

The treatment of the small scale physics of black 
hole mergers can be improved in several ways, for exam- 
ple considering the influe nce of the stellar density pro- 
file llVolonteri et a l.' '2003a") and gravitationa l kicks effects 
l|Tanaka fc Haiman .2009; ,Baker et al.1 l2008h which might 
cause the ejection of a fraction the BH from the host galaxy. 
Gas accretion could be modeled in a somewhat more real- 
istic way through the use of Bondi's formula; however, the 
validity of such approach is unclear. In addition it would 
require specific assumptions on the gas density profile. It 
would also be interesting to attack some of the unresolved 
aspects related to disk fragmentation induced by the pres- 
ence of metals. Finally, we neglected the radiative fe e dback 
of the quasars on the surrounding gas: IValiante et al.1 (|201ll ') 
showed that by means of this mechanism star formation in 
the host galaxy can be efficiently quenched and the gas ex- 
pelled. This process should become important only in the fi- 
nal phases of the evolution, when the accretion rates become 
extremely large and therefore it should not have dramatic 
effects on the high-z evolution with which we are mainly 
concerned here. 

Future observational developments that could in prin- 
ciple put some constraints on our model include the activi- 
ties of the International X-ray Observatory (IXO) for which 
we aim at making detailed predictions in the future, and 
those of the Laser Interferometer Space Antenna (LISA) 
that should be capable to detect gravitational wave bursts 
originating from bl ack hole merger events at high redshifts 
l|Sesana et al.ll20fll ). 
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